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We discover new vacuum solutions of the Jalilian-Marianlancu-McLerran-Weigert-Leonidov- 
Kovner equation, which correspond to center of a gauge group. We improve the color glass conden- 
sate (CGC) model by an explicit usage of a density matrix. Studying scattering of CGC states in 
an external color field, we observe that an amplitude is naturally expressed via group characters. 
We construct an example that shows how new thin effects may be potentially observed in peripheral 
collisions. We prove that at any parton density a gluonic CGC state does not become a true black 
disk. We find a wave function of a true black disk and show that it necessarily contains many 
quarks. This result corresponds to the necessity of nonvacuum Reggeon loops in a formation of a 
"q^ ■ true black disk. 

PACS numbers: 12.38.-t, 25.75.-q 
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t— 5 . I. INTRODUCTION 

o : 

Consider high energy processes in the QCD. A hadron at a given rapidity is a composite state of quarks and gluons. 
It is described by a wave function which is some state in the Fock space of partons. If we boost a hadron to high 
rapidity, then its wave function will be changed according to the methodology of evolution equation. When studying 
CIh' realistic processes, we usually work with an effective coarse-grained theory with a given scale of transverse resolution. 
JL , If we change the scale then a wave function will be changed, too. Then we can obtain the DGLAP equation (or usual 
■ RG flow for the Lagrangian, if we take a spatially symmetric case). Usually, the scale of transverse resolution is given 
, by a characteristic momentum transfer Q. Taking a more coarse scale, we shall calculate a wrong result. Taking a 
l— ~~ '■ more fine scale, we shall make unnecessary work because in calculations of observables we effectively average over 
| modes living far away from the interested momentum transfer. 

If at some transverse scale the number of partons per an elementary resolution area is much larger than 1, then 
we can expect an appearance of some sort of a statistically effective description. Therefore, the so-called color glass 
I/") , condensate (CGC) model [l[ was proposed, where specific conditions on operator correlators of color charge density 
was imposed. A next natural question is about an explicit partonic structure of the CGC model. We want to know 
what kinds of quantum states produce the CGC description for field correlators. There is a partonic derivation of 
\Q the CGC model for a large nucleus, or for an arbitrary dense state H. In the work fH, a large number of random 
color partons in the same place were considered. In our framework, we associate such a place with an elementary 
cell in the transverse plane. Reference Q concentrates on a distribution of a Casimir operator and its main goal is 
to justify the validity of the classical description, which exploits a classical color charge and a Gaussian distribution 
as a first approximation. However, such a description is too peculiar and is devoted only to the construction of an 
effective theory In this paper, we look at the situation from a more general viewpoint where a fast hadron projectile 
is a quantum system of QCD partons and a wave function is an element of the Fock space. We shall show that Ref. 
Q implicitly use a density matrix for a wave function of partons. In this paper, we emphasize the usage of a density 
matrix. This allows us to see properties, which lie beyond the scope of classical charges. In particular, we show that 
the CGC states do not tend to a black disk at arbitrary high parton densities. 

An important ingredient of high energy QCD is the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov- 
Kovner(JIMWLK) equation [J 0]. It describes the evolution of a scattering amplitude of a dilute projectile on a 
dense target. The problem of white and black discs was discussed in context of the JIMWLK/KLWMIJ equation 
It was shown that there are only two stable points of the evolution: black and white. In this paper, we show that the 
situation is more complex. On a compact group, instead of a usage of common Fourier transform we must use matrix 
elements of irreducible representations as a complete basic in a space of functions. This gives new results which were 
missed in Q. Particularly, we show that there exist multiple stable points of the evolution. In fact, they are the 
central elements of the gauge group. 

In non-Abelian gauge theories, the nontrivial role of the SU(N) center is already known. In the thermal QCD the 
spontaneous broken center symmetry is a property of the deconfinement phase transition (see Q for a nice review 
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and references therein) . The so-called Polyakov loop is an order parameter for center symmetry breaking. The group 
center, group characters, and the group manifold also emerge in various approaches to the problem of the confinement 
in QCD. For an example, in Ref. pj a random walk on the group manifold was observed. 

Our new results raise an important question about the structure of a true black disk. This disk is just a mathematical 
idealization. The precise definition will be given is Sec. IVIIII In Ref. [5j, a black disk was considered as a special 
choice of a weight functional which determines the probability density to find a given configuration of classical color 
charges in a projectile. The main flaw is the usage of classical charges. In this paper, we argue that classical charges 
are not suitable for question about a black disk. Our calculation contains only quantum operators which act in the 
natural projectile Hilbert space. We show that in order to construct the true black disk we must take many partons 
in the fundamental representation. 

The master plan of our paper is (1) To explore a partonic structure of the CGC model. (2) By symmetry reasons, 
to consider the CGC state as a projectile. (3) To study a scattering of the CGC in an external color field. (4) To 
investigate the limit of a high parton number. To find points on the group manifold where S-matrix does not vanish. 
(5) To try to find states which in high density limit tend to the true black disk. (6) To argue that a true black disk 
must contains many quarks. 

This paper is organized as follows. In Sec. UH we review basic properties of the CGC model and observe that there 
is a natural description of CGC in terms of a density matrix with maximal entropy. Then we show that an average 
of a gauge transformation on the CGC state is equal to the high power of a SU (N c ) group character. In Sec. IIII1 we 
argue why it is necessary to study a scattering of any state in an external color field. It will be clear why the method 
of weight functionals is wrong in a general situation. In Sec. IIVI we analyze a simple situation of the SU (2) gauge 
group and find points on a group manifold where the character of the adjoint representation has maximal value. In 
Section [V] we study a more complicated case of the ST/ (3) gauge group. In Sec. EH we show a scattering process 
where a center of a gauge group influences a physical observable. In Sec. IVII1 we prove that central points are stable 
solutions of the JIMWLK equation. In Sec. IVIIII we discuss a question about structure of the true black disk. Section 
IIXI contains our conclusions. 



II. CHARACTERS IN CGC 



Consider a hadron state \*f/} having a large number of partons. Let us take a small area Si within hadron and 
consider a subsystem in S± as a quantum system. It was motivated in [8( that in such situations it is useful to use a 
density matrix for the subsystem Si. Moreover, we can expect that in a chaotic environment the entropy tends to a 
maximally allowed value. Maximum of the entropy suggests that all available states are equiprobable and the density 
matrix is proportional to the unit matrix. 

Consider N partons in color representation R. The space of states H of the system is a tensor product of the spaces 
of state of every parton. 

H = y R ®V R ®...®V R = Vj N (1) 

N 

where Vr is a vector space of representation R. A color charge operator is a generator of the SU(N C ) group 

p a = T£ <g> 1 <g> . . . <g> 1 + 1 <g> T£ <g> 1 . . . ® 1 + . . . (2) 

where TjJ is the generator in the representation R. 

Taking the unit matrix in the space H in ([1]), the normalized density matrix of maximal entropy is given by 

w = -JW ( 3 ) 

a R 

where d R is the dimension of representation R. Averages of operators, which act only on H, are calculated as 



(6) = Sp(wd) = -Ls P (d) (4) 



'R 

For example, we can calculate the average of the quadratic Casimir operator C 



(C) = (p a P a ) = ^-Sp(T R T R ) = NC R (5) 

UR 
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where we used the obvious fact that Sp(T%) — 0. In the subsequent calculation we usually omit representation index 
R. 

It should be noted that, in general, we must also consider some distribution on N . In the statistical limit of large 
N it usually has Poisson-like behavior. Here we consider a simplified case with fixed N which is equal to the average 
parton number in the cell. 

Now we shall show how to obtain the usual CGC description in terms of classical color charges with the Gaussian 
distribution. We want to calculate averages of arbitrary polynomials of p a . Using identity Sp(T^) = 0, formulas @ 
and ([J|, we directly evaluate 

(P a ) = (6) 
(pV) = ^-Sp(T a T b ) = ■^ l S ab = p 2 8 ab (7) 

(p a p b p c ) = ^-S P {T a T b T c ) (8) 

(p a p b p c p d ) = M^ii [Sp{T a T b )Sp{T c T d ) + Sp{T a T c )Sp{T b T d ) + Sp{T a T d )Sp{T b T c )] 
+ -^Sp(T a T b T c T d ) 

In the large N limit we have 

(p a p b p c p d ) — > (p a p b ) ( P c p d ) + (p a p c ) (p b P d ) + (p a P d ) ( P b p c ) (10) 

and similar formulas for an higher even power of p a . In (|10J). we can easily recognize the Wick's rules. Higher 
correlators are reduced to all possible pair contractions. If we temporarily forget about odd powers, then we can 
generate averaging rules with the help of the classical Gaussian distribution, which is described by the following 
weight functional: 

/ \ N *~ l 

^ l= bW " p (-2^' v ") < n) 

The average classical charge density we define as an average of the absolute value of the following classical variable: 

(Kl) = j \p a \W[p]dp ~ VN (12) 

Usually, the square root behavior is interpreted as the random walk scenario. The average of the gauge transformation 
is 

(e^ a )=exp(-ip 2 a Q a Q ) (13) 

For large N in (|13p . we have more and more narrow peaks near a a = 0. 

Now consider odd powers of p a . It is clear that there is a direct analog of (flQ)) for odd powers in the large TV limit. 
Namely, at odd n the leading term has only one triple contraction ©. 

(p ai ...p an ) = (p ai p a2 p a3 )(p a4 p^) . . . (p 11 ^ 1 p a ") + permutations (14) 

Note that (p ai p a2 p°" 3 } in general is neither symmetric nor antisymmetric. Hence, we cannot construct an analog of 
(1111) and cannot use classical variables. In order to cure the problem the so-called Wess-Zumino term can be used 
9, 10]. Here we use a different method of generation of averages based on matrix elements of gauge transformations. 

Now we switch to the new CGC description which uses functions on the group manifold. In Eq. {2j the operators 
p a are generators of gauge transformations. On the other side, the gauge group action on a tensor product is a tensor 
product of linear operators which represent group elements 



U A ® B (g) = U A {g)®U B {g) 



(15) 
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where Ua{9) is a matrix which represents g € SU(N C ) in the representation A. Since any matrix U(g) is a quantum 
operator in a color space, we can evaluate its average by applying the definition ^ with the density matrix ([3]). This 
gives 

%] = (U(g) ® U(g) . . . ® [7( 5 )) = " (16) 

where Xi?(g) is a character of the representation R. Characters can be considered as invariant functions on a group 
manifold. Basic properties of characters will be reproduced in Sec. IIXI For SU(N C ), group characters are uniquely 
defined on diagonal matrixes. Any diagonal matrix is defined by its eigenvalues z%. Unitarity leads to \z,i \ = 1. Hence, 
we have constraint 

xn(g) 



< 1 (17) 

So it is clear from (fTE)) that in the large N limit only the neighborhood of points \xr\ ~ <2r have a nonvanishing S 
matrix. Moreover, if \R = dn e where 0^0, then due to existence of a distribution over N, as it was discussed 
above, such points do not survive in the large N limit. So we are primarily interested only in points where \R = 
and where S ~ 1 , even in N — > oo . 

Let us introduce canonical coordinates a a on the group manifold. Below we shall simultaneously use two notations 
for group points: a set of coordinates a a for an analytical notation, and an element of group g for an invariant 
notation. Using the coordinates, an average of a gauge transformation has the form 

S[a] = (e ia «» a ) (18) 

The group multiplication induces two kinds of global vector fields on a group manifold, the so-called left and right 
invariant vector fields: and J". They act on S as 

JlS[a] = i(p a e iabpb ) (19) 



J a _S[a] = i{e la " pb p a ) (20) 

An average of arbitrary power of p a can be calculated by the action of J± on function S[a] at the point g — e (or 
a a = 0), where e is a identity element of a group. 

Physical meaning of S[a] is an elastic S matrix of a scattering of a CGC state in an external color field a a . Consider 
a collision of the CGC state considered as a projectile with some other state named as target. The result (|13l) is valid 
only near a a = because (fTU)) is valid only if power of p a is much less than N. Moreover, (fT3")l is not a function on 
the group manifold. Our gauge group is compact and has finite volume. Any function must be at least periodic in 
the context of canonical coordinates. The correct value of S[a] for CGC states was calculated in (IT6l) . In addition, 
the result (|16[) is valid for all values of N. Relatively large N is required only for a search of regions where S[a] is 
maximal. 

Finally, one more note is about the definition of CGC. Originally, it was destined to describe dense states with a 
large number of partons per an elemental cell. However, it is clear that the CGC model is not rigidly restricted to 
the dense domain. It is only required that all correlators must vanish, instead of 

(p a (x)p\y)) =p 2 5(x-y) (21) 

Partons at different transverse cells must be decorrelated. Using a density matrix, it is clear that such decorrelation 
naturally arises in states with large or maximal entropy. If we select some subsystem, then we can use entanglement 
entropy. A small dilute subsystem, as a part of a large dilute system, may be in a state of maximal entropy due to the 
combinatoric reason. There is also a question about the case x — y in (f2"Tj) . We show now that the formula (|2"Tj) may 
be rewritten in a more natural form, which has no any subtleties like x = y. The right-hand side of (|21[) is entirely 
related to yields of the normal ordering procedure. Let us recall that 

P a = r'u.u, (22) 

We wish to know the operator product p a p b , where the operators act on same transverse position. Performing normal 
ordering, for the fundamental representation we obtain 



p"p>> = : p y : +C abc : P c : +^ ab N 



(23) 
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where C a bc is an unessential group tensor and N = ^2 a \ a i 1S the operator of parton number. The average (N) gives 
corresponding /i 2 factor in (|2"Tj) . Hence, the CGC model can be reformulated as the following requirement for the 
state |*): 

(*| :p ai ...p a " : |*) =0 (24) 
for any n > 0. The formula (|24|) is a natural replacement for (|2T1) . 



III. GENERAL PICTURE OF HIGH ENERGY SCATTERING 



In order to motivate the necessity of studying scattering in an external color field we propose a general picture 
of high energy scattering. Consider an arbitrary projectile-target scattering. We use Cartesian coordinates with the 
instant form of dynamic, which allow us to write an explicitly symmetric theory. The interaction occurs only in a 
short time interval near t = 0. Let |P) be a state of the projectile and |T) be a state of the target. The initial state 
at t — — oo is 

|P> \T) (25) 

Note that this state is a pure tensor product. We want to know the amplitude to find the system at t = oo in the 
given final state |P') ® \T') which is also a pure tensor product. The amplitude is a matrix element of the quasielastic 
S matrix and it is given by 

S = (P'| <8> <T'| S \P)®\T) (26) 

where S is the scattering operator. Any initial state is a quantum superposition 

|P)=$>f|P fc ) (27) 
fc 



\T)=J2*l\Tk) (28) 

k 

where \Pk) and \Tk) are states of pure quasiclassical fields without a quantum interference. Also, they play a role of 
boundary field configurations for path integral at t — ±oo. In a theory of point particle, such states are just the set 
of wave functions 5{q — qo). In a field theory such states have a fixed field configuration and a wave functional has 
a form 5{ip — ipa). In the high energy QCD, as a source of bremsstrahlung field a fast parton with a fixed transverse 
coordinate generates a classical gauge field without interference. Vice versa, this parton propagates in the target field 
considered as an external field for the parton motion. If the parton has a smeared wave function, then according to 
the superposition principle its total amplitude is just a linear combination of the quasiclassical amplitudes which is 
usually expressed via Wilson lines. Note that an interference in color indexes can be encoded in terms of classical 
configurations because the classical field also has a color index. However, the high energy evolution converts any 
classical state to a state with an interference, since new emitted gluons arise with a spreaded wave function. A state 
described by a density matrix also does not generate a pure field. 

Following to the usual quasiclassical ideology, we express the action of S on any initial state via a superposition of 
same action on a pure fixed quasiclassical state. In the quasiclassical approximation an amplitude is just the exponent 
e lSa where Sq is value of a classical action evaluated on a classical solution of motion equations which give extremum 
conditions of the classical action. So the action of the S matrix has a form 

S \P)®\T)=J2*k^S P (aT)\P k )®S T (a?)\T s ) (29) 

fc,s 

The projectile propagates in the target fields and the target propagates in the projectile fields. Equation (f2"9"| is valid 
since in a high energy scattering the time interval of the interaction is very short and a change of the target state does 
not affect on target fields viewed from the projectile viewpoint. After the collision, causality disallows an influence of 
the changed target on the projectile. The field a? is a classical field generated by the target state |T S ). In the QCD 
it is a gauge field in the Lorenz gauge emitted from classical parton currents . Note that the final state in (j29|) 
is not a tensor product. It cannot be expressed as |P) <g> |T). 
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To find an elastic amplitude we set \P') = \P) and \T') — \T). It has a form 

S el = ]T *fa^^^( P fe|5p(<)|P fcl )(T S2 |,STK)|T S2 ) (30) 

ki ,k2,Si ,S2 

If we take an asymmetric case with a dense target and a dilute projectile, then is small and Eq. (|30|) converts to 

S el =J2(P\Sp(^)\P)\^\ 2 (31) 

s 

This is the well-known representation for an elastic amplitude. The projectile scatters in the external field and the 
result is averaged over various target field configurations with the weight given by target wave function For 
a general case, Eq. (|3U|) cannot be expressed in the form of Eq. (|3"Tj) . The conception of weight functionals is 
hardly inapplicable because the final state is the entanglement state which is not a tensor product. Nevertheless, 
the amplitude (Pk 2 \S(a T )\Pk 1 ) is the main building block for the exact expression ([50)1 . This amplitude for fixed 
projectile states gives a functional over target field. 

Concerning the question about a black disk, we define a black disk as a state \B) which for any nonempty projectile 
\P) gives exactly zero elastic amplitude. The usual way is to apply the assumption that \B) is dense and \P) is dilute. 
In this case, the black disk condition is converted to a condition on the weight functional | 2 . Indeed, the standard 
CGC model at very large square charge density poses a good approximation to the black disk. However, the restriction 
into the dense-dilute case is too crude. The true black disk is defined by the more strong condition. Namely, S e i in 
Eq. (|30|) must be equal to zero for any nonempty state \P). Note that this requirement is projectile-target symmetric 
by the construction. Currently, we are unable to treat directly the general case. So we restrict ourself to a case where 
\B) and \P) are classical states which were defined previously. In this case, the black disk condition is 

(P\S P (a B )\P) (B\S B (a P )\B) =0 (32) 

which is explicitly symmetric, again. 

To satisfy the requirement (|3"2"|) we cannot claim (P\Sp(aB)\P)—0 for any \P) because such ap does not exist. Only 
the second condition can be satisfied 

(B\S B (a P )\B) = (33) 

for any ap ^ 0. This is exactly why we have studied the scattering of the CGC states in external color fields. Further, 
by finding \B) which satisfies Eq. (f3"3")> . in Sec. I Villi we shall explore the structure of the true black disk. 



IV. SU(2) GAUGE GROUP 

Consider the SU{2) gauge group as a simple example. The diagonal subgroup has one dimension and is generated 
by the element T z 

°) (34) 



e-*2V 



where tp € [0;47r). In the adjoint representation we have 



The character is 



e t(p 

U 3 (g(<p)) =(01 ) (35) 
e- i( f 

Xz{f) = 1 + 2 cos ^ (36) 

Here we see that X3(2tt)/3 = 1. This point corresponds to — 1 € SU(2), which is an element of the center subgroup 
Z2. Hence, in some large nonzero external field corresponded to —1 element of group, a CGC state constructed from 
many SU(2) gluons can scatter without interaction at any gluon density. And moreover, in Sec. IVIIi we shall show 
that the state whose S[a) is localized near the point g = — 1 is a stable point of the JIMWLK evolution equation. 

In general, the irreducible representations of the SU{2) group is labeled by index I = 0, 4,1, | — Characters can 
be directly calculated as 

sin (I + |) <p 

Xl(<P) = . 1 (37) 

sin 

In general, the condition x;(27r)/d; = 1 holds at any integer I. 
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V. 517(3) GAUGE GROUP 

In this section, we shall show how to calculate characters of the SU(3) group. A first-time reader may skip this 
section without a loss of the paper completeness. 

The general situation can be considered with the help of the Cartan classification of Lie algebras. A character is an 
invariant function on a group: %(</) = xi^^gh) for all g,h € SU(N C ). For any g we can find an element h such that 
the element h~ 1 gh is diagonal. All diagonal elements form an Abelian subgroup, Lie algebra which is the equivalent 
of the Cartan algebra H. Often this group is called maximal torus. Irreducible representations can be classified by 
studying vectors in the dual vector space H* which is a linear space of linear forms on H. For each representation 
there is a special finite set of forms which are named weight vectors. For the adjoint representation (the Lie algebra 
itself) such weight vectors are named roots. In Fig. [TJ the fundamental and the adjoint SU(3) representations are 

\ /(§.!) 
?(°>f) \ 

(1,0) 




FIG. 1: Fundamental and adjoint 517(3) representations. 

shown. Bold points are the weight vectors. Roots can reproduce properties of all other irreducible representations. 
Let v be a vector in some representation and T(h)v = u(h)v, where w is a weight and h G H. If a is a root, then 
T(h)T{E a )v = T([h; E a ] + E a h)v = (a(h) + uj(h))T(E a )v. Hence, the vector T(E Q )v has weight a+wor zero. In 
addition, we can see that roots set length of all weights. Weights are associated with eigenvalues of T(h). So we can 
directly calculate the character on the diagonal subgroup 

XR (h)=J2^ kih) (38) 

k 

where the implicit exponential map from H to the diagonal group was assumed. For SU(3) in the fundamental 
representation we choose the following paramctrization: 

, ( -I 0\ 

(39) 

The element h\ corresponds to the x axis and hi corresponds to the y axis in Fig. 1. The diagonal subgroup is 
generated by 

g{ip, i>) = ex.p(i(phi + iiph 2 ) (40) 
where tp € [0;47r) and ip £ [0;67r). The character of the fundamental representation is 

X3 (v,V')=e i2 ' / 3 + 2cos|e-^ 3 (41) 

In the adjoint representation, we have eight vectors. Two vectors have zero weight. Six others are generated by 
symmetries from two positive vectors: (1; 0) and (^; 1). Character of the adjoint representation is 





2 + 2 cos ip + 2 cos (| + tp) + 2 cos (f - tp) 



(42) 



In a general case, our calculation is generalized to the well-known Weyl character formula. Now we want to find points 
where xs — 8. There are six pairs (ip,ip), where this condition is satisfied 

(0,0) (0,277) (0,4tt) , , 

(27T,7r) (27T,37T) (27T,57T) [ 0) 
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Substituting these points into (140]) , we obtain only three independent diagonal group elements: 1, e l27r / 3 , e l47T ^ 3 . Indeed, 
they form the center subgroup Z 3 . Geometrically, the map from a group to a manifold of linear operators can be 
viewed as n-fold covering. Points where \R = correspond to the kernel of the map. In the SU (2) case, we have 
the 2-fold covering SU(2) -> SO(3). 

VI. AN EXAMPLE OF TUNNELING 

In order to give a physical manifestation of the explored feature, in this section, we shall construct a gedanken 
scattering process whose elastic scattering amplitude is a black disk with "holes" corresponding to the central elements 
of the gauge group. For the sake of simplicity we consider the SU(2) gauge group. Let the projectile be a small (in 
the transverse plane) gluonic CGC state whose average parton density is not required to be high. As it was shown 
in Sec. IIII1 to calculate the full elastic scattering amplitude by using the formula (f5Tj) we must assume that the 
projectile field is small at points where the target is located. Simultaneously, to see tunneling the target fields at the 
projectile points must be arbitrary large. These requirements may be achieved by the following methods: choosing a 
dense large target, choosing a small dilute projectile, and increasing the transverse distance between the target and 
the projectile. One more simplifying assumption is that the target field oi^(x) has a distribution \^^\ concentrated 
strictly near only a one classical field configuration. In this case, the formula (|31[) has a simple from 

S(x) = (P\S P (a T )\P) (44) 

where in the evaluation of the matrix element we can neglect x dependence of a T if a function a T (x) is slow varying 
at the projectile points. Since \P) is a CGC state, we have the picture discussed in Sec. [TT1 where a CGC projectile 
scatters in a constant external color field. Hence, for SU (2) gluons we can guess 

S(x) = (jM*)))" (45) 

where N is a number of gluon in the projectile, g(x) is the group element corresponded to the group coordinates 
a I(x)- 

Strictly speaking, Eq. (|45[) gives a nontrivial contribution only if the projectile is not an overall singlet state. 
Obviously, a singlet state in a constant external field has exactly S = 1. Nevertheless, there are no reasons to drop 
nonsinglet states because a general formalism of scattering must be able to treat the whole set of states on same 
mathematical ground. In principle, to make the projectile singlet we can add at large x a small patron cluster that 
discharges the projectile. Since at large x the color field is small, new partons have 5=1 and do not change the 
projectile S matrix. The situation can be shown in a more straightforward way. Let the projectile be a large gluonic 
dipolc. Its S matrix is 

S(y,x,a) = -x 3 (g~ 1 (y)g(x)) (46) 

Let x be fixed. Let us move y to a point where a ~ or to the infinity. This allows us to approximate (|46|) as 
S ~ X3(g(x))/3- So, dealing with a singlet state, we again obtain a expression similar to (|4"5| . Physically, this is a 
special kind of a border collision where one dipole's parton probes the target and other is a outlying spectator. 




X 



FIG. 2: Left panel: The scattering setup in the transverse plane. The projectile is a small gluonic CGC state. The target is 
a dense composition of localized patrons. The target field a a (x) is fixed. The impact parameter x is large. Right panel: The 
schematic plot of S(x). The narrow white peaks arise in the black area. 



We discussed in Sec. IHII that a parton with a fixed transverse coordinate and a fixed color state generates a classical 
gauge field without an interference. Two such partons generate a pure field, too. So we construct the target as a 
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dense color multipole. Namely, the target is an aggregation of a large number of independent partons which generate 
large, undistributed total field at the projectile. We collect our scattering setup in the left part of Fig. [51 

To calculate the scattering amplitude (|45|) we need to know the target field aj(x). We estimate it as a color 
multipole which at large distances has a power-like asymptotic. Taking a fixed transverse direction of the vector x, 
we guess 

where x = \x\, and c a is a color factor which describes the color charges in the target. By unitary transformation in 
SU(2), the matrix c a T a can be diagonalized as cc 3 /2. Hence we can apply Eq. (f3"B| which gives 

S(.)-(l±2=^)" (48) 

A schematic plot of S(x) has been drawn on the right part of Fig. O At large x there is a white area S = 1. However, 
at low x (black area) there are narrow white peaks which are a new result. The original CGC model has no such 
peaks. The origin of these peaks is the central elements of the gauge group. When the field a T is near a central 
element, the scattering amplitude fast rises. A physical consequence of this result is that at some exceptional impact 
parameters the inelastic processes is strongly suppressed, since there is only elastic channel at S = 1. This explains 
the exploited term "tunneling" . The phenomenon occurs only for a non-Abelian gauge group. For an Abelian case, 
characters always obey \\\ = 1. 

VII. NEW VACUUM SOLUTIONS OF JIMWLK EQUATION 

The JIMWLK evolution equation 0] can be expressed in terms of invariant vector fields on a group manifold [TT[ | . 
We also use the notation taken from Ref. [ll|. The JIMWLK equation is 

^ = J K xyz {-r + (x)Jl(y) r{x)r_{y) + 2V ba {z)j\{x)J a _{y)) S[a] (49) 

zxy 



9 2 jz-y)(z-x) 
(2tt) 3 {z-y) 2 (z-xf 



where Vb a (z) is a gluon scattering amplitude in an external field a a . We have shown in the previous section that at 
the points go : X-r(So) = <^R the operator Un(g) = 1. Consider a model functional 

S[a}=l[5(g(x)-g ) (51) 

X 

where the expression g — go is formal and has sense only in coordinates. If S[a] is localized near go in the whole 
transverse plain, then in Eq. (|4"9"| we have Vb a (g) — > Vb a (go) = $ba- The left and right invariant vector fields are 
related to each other by the adjoint representation 



Jim = ^-Jihg) 



= ^f{gg- l hg) 



V a b(g)J b -f{g) (52) 



h=0 



dh° 

Using the fact Vb a (go) = $ ba once more, we have 

JtS(g - go) = rS(g - g ) (53) 

where instead of S function we can take any function localized near go- Substituting property (|53[) into (|49p . we 
see that S(g — go) is a solution of the JIMWLK equation with zero eigenvalue. Strictly speaking, such solution is 
not physical because a physical state must obey condition S[0] = 1. Though, physical states may have a nonzero 
projection on S(g — g ). For example, in Sees. IIVI and fVl we have seen that the gluonic CGC state indeed has the 
nonzero projection. In relation with Ref. [f|, where an existence of only two physical vacuum states: black and white 
was motivated, our new states are grey due to a possibility of tunneling through some external nonzero field. 
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VIII. TRUE BLACK DISKS 

In the previous sections, we have shown that the CGC model does not contain a black disk, since in the special 
external fields any CGC state has a nonzero S matrix. Next, we shall step beyond the CGC model. It is natural to 
ask what projectile states have an exactly zero S matrix on a whole group manifold instead of a unit element. Such 
state is a natural candidate to a black disk. 

Usually, a black disk is defined as some dense state which absorbs any projectile. In this section, we consider a 
reverse case where a black disk is a fast projectile and scatters in an external color field. Any hardon in the QCD 
can be described by a wave function represented as certain vector in the Hilbert space. Consider a collision of a black 
disk state |^) with a some target. We take a small black disk and a very large and dense target. We suppose that the 
target color field is much stronger than the projectile field. This does not suppose that the target is a black disk. The 
key issue here is that, freely choosing target, we can construct an external color field as we want. In the considered 
case, the projectile elastic S matrix has the form 

S = J (V\e l f aa{x)pa{x)dx \y)W[a]Da (54) 

where a a (x) is a classical field generated by some quasiclassical element of the target wave function, p a (x) is the 
operator of color charge (the generator of gauge transformations) which acts on the projectile state, x is a transverse 
position. The functional W[a] is a weight functional which is obtained from the target wave function. Since the 
target is arbitrary, we shall not perform the average on a a and consider a configuration with fixed a a . So the elastic 
projectile S matrix equals to the functional S[a]. 

In the previous speculations, there was one mistake. Though the field a a {x) can be very large, due to compactness 
of the gauge group the scattering amplitude is bounded. So, in general, we have no way to compare fields a a (x) in 
the context of scattering. Nevertheless, we define a black disk state as a state which has a zero S matrix if a a (x) 7^ 
in the transverse area where the black disk is located, and which has S — 1 if a a {x) — in the same area. We shall 
name this state as true (or extremal) black disk because usually in current literature more weak definitions are used. 
We define 

S[a] = {^\e$ a " {x)pa[x)dx \^) (55) 

where it should be stated that the x integral in (|55[) is bounded by points where p a {x)\^>) ^ 0. Since the true black 
disk must be black for any target and any Wfa], we must claim 

S[a] = z\{5{a a {x)) (56) 

x,a 

in the sense of functional generalization of the 5 function. The factor Z is the normalization constant which is not 
important for us. It was proposed in [j| to use variables u) a (x) which arise from Fourier transform of S[a] 

S[uj] = J S[a}e l I aa{x)uja{x)dx Da (57) 

For the infinite black disk this gives a constant functional. The variables w a (x) are usually called classical color charge. 
But this method is applicable only for small a a . For arbitrary a a we should not forget about SU(N C ) compactness. 
Via the exponential map e la " p (a map from a Lie algebra to a representation) the variables a a (x) can be viewed as 
canonical coordinates on a group manifold. The integration over a a {x) must be considered as the integration over 
a group manifold over Haar measure 1 . Instead of a a {x) we introduce the variables g(x) which is a map from the 
transverse plane to the group manifold g(x) : R 2 — ► G. The claim (|56l) can be rewritten as 

S[g}=Y[6(g(x)-e) (58) 

X 

where e is the identity element of the group G. The function 5{g — e) can be viewed as a function on the group 
manifold with localized suprcmum. The expression g — e is rather formal and has a sense only in the coordinates. In 
general, we define 8 function as a linear functional over functions on a group manifold: J S(g — e)f(g)dg — /(e) 



1 For the diagonal Killing form Sp(T a T b ) ~ 8 ab and small a a the Haar measure degenerates to the usual form Y\ a dot a - 
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On a compact group instead of using the Fourier transform we must use the matrix elements of the irreducible 
representations as a complete basic for space of functions. The well-known Peter- Weyl theorem says that function 

b%(g) = v^yG?) (59) 
forms an orthonormal and complete basic in the Hilbert space L 2 (G, dg): 



bmKj(g)dg = SrqSmS^ (60) 

The letters R,Q denote irreducible representations; dn is dimension of a representation; and 7r^(<?) is a matrix element 
of the representation R on the group element g. 

A function /(<?) is called central function (or class function) if /(hgh^ 1 ) = f(g) for any h E G. It is easy to show 
that the function 8{g — e) is central due to invariance of the Haar measure 

Sihgh- 1 - e)f(g)dg = f S(g - e)f(h~ 1 gh)dg = /(e) (61) 



A direct consequence of the Peter- Weyl theorem is that characters of irreducible representations form an orthonormal 
and complete basic in the Hilbert space of central functions. Characters are given by 



Xr(9)=J2^(9) (62) 



f XR(g)XQ(g)dg = 6 R q (63) 

Any central function /(<?) can be decomposed in series on characters 

/(s) = X^«XKte) (64) 

R 



fn = J XR{g)f(g)dg (65) 

Here we are interested in a decomposition of the 5 function. 

^9-^)=T, d ^R(g) (66) 

R 

where we used the obvious fact that Xi?( e ) = dn- It should be emphasized that in a case of the gauge group SU(N C ) 
at N c > 2 there are representations which are not equivalent to its complex conjugate. For example, a quark is not 
an equivalent representation for an antiquark. Such representations have complex characters, but they are not purely 
imaginary near e. The formula (|66[) is a rather mathematical idealization. Infinite series are irrelevant for the physical 
reasons. In the real world, the S function must be smeared near e. So the black disk S matrix is not zero for very 
small fields a a . This is the equivalent of a truncation of the series (|66p which, at first sight, is hardly divergent due 
to the dji growing. The convergence of (|66p is achieved by alternating characters values at nonsingular points. High 
representations in ([6"6]) correspond to high frequencies in the usual Fourier analysis. This happens due to the growing 
of the second Casimir operator Di when dimension of a representation grows Q. The second Casimir operator on 
functions is the Laplace operator on the group manifold. The eigenvalues of the Laplace operator correspond to 
square of frequencies vector in a flat case. It well known from the wave packet calculations that this high frequencies 
suppression is equivalent to spreading of a wave packet. In order to obtain a more narrow packet we must include 
higher harmonics. As an example, in Sec. IIV1 the characters of the group SU (2) were calculated explicitly. 

Characters have useful properties. Cyclicity of the trace in (|62[) gives x{hg) = x(gh)- The group multiplication 
induces two kinds of global vector fields, which are called left and right invariant vector fields: J" and J*. Namely, 
there is a linear map from Lie algebra to vector fields on the group manifold. Acting on characters, the invariant 
vector fields give equal results. 

= J-x(g) (67) 

h=0 



J a +x{g) = gj^xihg) 



- O^x(gh) 
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From unitarity of representations we have 



J a ±x(g) 



= o 



(68) 



In Sec. IVII1 we have shown that the functional ([56]) obeys the JIMWLK equation. This means dS[a]/dY = 0. 
Strictly speaking, the JIMWLK equation is valid only in the dilute regime, so our calculation has only demonstrative 
purposes. 

The JIMWLK Hamiltonian is a linear operator. If we take its eigenvalue —A, then at high Y the corresponding 
eigenfunctional disappears as exp(— \Y). The true black disk has zero eigenvalue. So it is a stable state of the 
evolution. This is not true at the black disk border where there is a transition region (crossover) to the white state. 
Obviously, any black disk must grow when we go to high Y. Hence, at the black disk border an interesting process 
will occur. Emitted quarks and gluons will saturate the border by higher characters, thus the border will become 
black. 

The next question is about a possible partonic interpretation of the characters xr- The most obvious method is to 
use a density matrix, again. Consider one parton of color charge R. Let the density matrix of the parton color space 
be proportional to 1 (maximum entropy state). The normalization factor is simple . Then the elastic scattering 
amplitude in an external field a a is 



where the averaging is assumed over the density matrix. The average color charge within the black disk is zero: 



This behavior corresponds to the behavior of CGC states, where, in addition, it is required that (p a (x)p b (y)) = 
for i / i/. In difference from the CGC model, in the our black disk state the averaged charge square p a (x)p a (x) is 
divergent. 

It follows from (|66[) that it is absolutely necessary to have many quarks and antiquarks in the black disk wave 
function. Many representations cannot be obtained via a reduction of tensor products of the gluon representation. 
For example, in the SU(3) case such representations are 3, 3, 6 . . . Currently, we are not able to predict how such a 
scenario actually can be realized in nature, but the necessity of existence of partons in the fundamental representation 
in the true black disk state was proven. So we can suggest the following hypothetical scenario. Starting from a dilute 
projectile such as a dipole, when boosting it to high Y we predominantly generate many new gluons in the wave 
function. When the gluon density becomes high, then there will be a large probability of emitting soft quarks into the 
new opened phase space. This corresponds to the mechanism of nonvacuum Reggeon in QCD. See [1 31 ] for a modern 
example of a quark propagation to a low rapidity. Usually, such processes are kinematically suppressed in comparison 
with a gluon emission, but we can expect that in a case of a high number of sources it can be very probable. A fast 
gluon converts into a fast quark and a slow antiquark and vice versa. 

Another possible interpretation comes from the running coupling effect studying (l4| . At NLO there are diagrams 
with quark bubbles. Usually, such diagrams are associated with change in the effective QCD coupling constant. 
However, in the context of the current paper, we see that such diagrams can really saturate the black disk wave 
function by the required fundamental representations. So the considered true black disk really may be an endpoint 
of an evolution equation constructed with higher order contributions of the perturbation theory. 

The fact that in order to obtain appropriate characters we took the density matrix of maximum entropy naturally 
corresponds to an expected chaotic behavior of a wave function boosted to high energy. 



In this paper, we have found the new vacuum solutions of the JIMWLK equation and have explored the structure 
of the true black disk. Although it is hard to find any relation of this result to available experimental observables, it 
gives a new theoretical insight to the QCD asymptotic at very high energies and densities. In perspective, we can use 
geometry on a group manifold as a powerful mathematical tool, which helps to improve our understanding of high 
energy QCD. In particular, we have shown that the CGC states do not tend to a black disk. We have shown that 
the true black disk has a more complicated structure and it must be constructed from the whole set of irreducible 
representations which may come from a composite state of initial partons. 

One can be confused by the fact of existence of a large number of quarks in the black disk wave function. However, 
we stress here that this notion is about an ideal black disk, which is defined by the rigorous mathematical definition 




(69) 



(p a ) = J a 5(g-e)\ g=e = 



(70) 



IX. CONCLUSION AND DISCUSSION 
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([56]) . Such black disks may arise only at infinity high Y as a final state of the evolution. Far before infinity Y we 
shall deal with some approximations, which are experimentally indistinctive from the true black disk. 

Let us discuss a relation of our method to a saturation scale. In a description of an experimental data a so-called 
saturation scale Q S (Y) is widely used. It is a useful quantity that delimits the high and low density regions of high 
energy QCD. Experimentally it is related to a dipole scattering amplitude. A dipole size, at which an amplitude 
becomes saturated, has the scale 1/Q S . Before the saturation there is a geometrical scaling where some observable 
depends only on the dimensionless quantity Q 2 /Q 2 S (Y). In the original CGC model Q s is associated with the average 
square charge density [i 2 , which enters in (|T|), (|13[) . Because of the evolution, the parton density grows and Q s grows, 
too. In the context of the current paper we used the clear partonic interpretation of a saturation scale. Namely, this 
is a scale, at which the average parton density is near 1. Quantitatively, it is proportional to square root from the 
transverse parton density, which is taken at some initial scale Qo (which is needed to reflect the DGLAP evolution). 
If we go to high Y , then the parton density will grow and Q s will grow, too. Making our conclusion, we see that Q s 
is too crude a characteristic of processes. To see more thin effects we have to use more susceptible quantities. 
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